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The concept of critical points in nuclear phase transitional regions is discussed from 
the standpoints of Q-invariants, simple observables and wave function entropy. It 
is shown that these critical points very closely coincide with the turning points 
of the discussed quantities, establishing the singular character of these points in 
nuclear phase transition regions between vibrational and rotational nuclei, with a 
finite number of particles. 
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Nuclear structural evolution in transitional regions is often thought of as a 
continuous variation of properties, as a function of nucleon number, from one 
idealized limit (e.g., vibrator, rotor) to another. The rapidity of structural 
change may vary across a transitional sequence of nuclei, and different mass 
regions exhibit different rates of change but, until recently, no individual point 
along these evolutionary trajectories could be singled out with special obser- 
vational properties. 

In the last years, however, the concept of critical points in shape/phase tran- 
sition regions has been much discussed [1-5]. While the concept itself is well 
known in nuclei (in the context of the coherent state formalism [6,7] of the 
IBA model [8]), it is only very recently that analytic descriptions of critical 
point nuclei have been given [9,10]. This is a significant point since, histori- 
cally, such nuclei have been the most difficult to treat: they exhibit competing 
degrees of freedom, and one has had to resort to numerical calculations. 

Two critical point symmetries, called E(5) and X(5), have been proposed 
[9,10], giving analytic expressions for observables which are exactly at the 
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critical points of a vibrator to axially asymmetric (7-soft) rotor transition 
region, and of a vibrator to symmetric rotor transition region, respectively, 
for an infinite number of nucleons. An important aspect of this is that, for 
the first time, one is able to associate special observational characteristics to 
a specific point along a trajectory from one structural limit to another. Re- 
cently [11], using the methods presented here, the well-known 0(6) limit of 
the IBA has also been identified as another, heretofore unrecognized, critical 
point symmetry, for the transition between prolate and oblate nuclei. This is 
an important result since the 0(6) symmetry can be calculated in the IBA for 
finite nucleon numbers, in contrast to the non-IBA symmetries E(5) and X(5). 
So far only two examples for nuclei [12,13] which lie close to the X(5) and E(5) 
symmetries are known while, interestingly, there are many examples for 0(6) 
like nuclei. In the present work we will restrict our discussion to prolate nuclei. 

To understand the evolution of structure in real nuclei, with a finite number 
of nucleons, it is important to gather information about systematic changes 
of observables at or near such critical points. This aim can be achieved by 
the use of a model that is able to describe limiting cases of nuclear structure 
- vibrators, rotors and 7-soft nuclei - and a large variety of nuclei between 
these limits. Such a model is given by the IBA, which - in the expansion of 
the coherent state formalism - exhibits critical points as has been discussed in 
refs. [6,7,14,15]. We stress that the critical point descriptions X(5) and E(5) 
are defined in terms of a geometrical approach, not the IBA. Nevertheless the 
IBA provides a convenient tool to span a range of structure, including phase 
transitions, and also to assess effects of finite particle numbers. 

It is the purpose of this Letter to show, from several complementary theoretical 
approaches, that there is independent evidence for the singular character of 
these critical points, and independent ways of identifying them in observables 
calculated in collective models. To do so we bring together three major themes: 
the already mentioned study of phase transitional regions and critical point 
nuclei, the behavior of quadrupole shape (Q)-invariants, and the study of chaos 
and entropy in nuclear systems. We show that the critical points occur very 
near to the turning points (points of steepest descent or ascent) of these Q- 
invariants - that is, at the extrema of their first derivatives. The same behavior 
will also be shown to hold for some more easily accessible observables. 

To span the transition regions, it is convenient to use the IBA Hamiltonian in 
the following form 

H = a[(l-()n d -^Q.Q} (1) 

where Q = s^d + d^s + x\d d}^ and we consider the well known parameter 
space of the extended consistent Q formalism (ECQF) varying ( between 
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and 1, and x from to —\flj1 = —1.32, while a is a scaling factor. This 
parametrization is equivalent to the more commonly encountered (equivalent) 
ECQF [16,17] form of H, which includes the parameters e and k. 

Figure 1 illustrates the three dynamical symmetries of the IBA in terms of a 
triangle. With the Hamiltonian of Eq. (1) it is easy to calculate the structure 
for any point in the triangle. For ( = one obtains a U(5) structure (for any 
X), and ( = l,x = — a/7/2 gives SU(3). Thus, a U(5)^SU(3) transition region 
is defined by x — /2 and ( varying from to 1, while a U(5)^->0(6) region 
has x — and ( varying from to 1. 

One can use the coherent state formalism [6,7] of the IBA model to identify 
the critical points in the ECQF space. In this approach, the energy functional 
for the ECQF Hamiltonian is given by 



*(C X, I 3 ' 7) = 1 + 02 * >~ . 2 . 



The variation of £ changes the structure between the vibrator limit and ro- 
tational nuclei - both axially symmetric and axially asymmetric - which are 
the transitions we will focus on. Critical points in ( are found where E be- 
comes flat at P — 0. These points, which we refer to as ( c , can be derived by 
evaluating the condition 



d 2 E{Q 
dp 2 



= 



(3) 
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On the transition path from U(5) to 0(6) (for x — 0) exactly one critical 
point is found, namely where a second, deformed, minimum in P of the energy 
functional emerges. 

The situation becomes more complicated for transitions with % 7^ 0. In these 
cases, the spherical minimum is joined by a deformed minimum and both 
minima coexist in a very close parameter range in (, converging to one point 
when approaching x = 0. Thus, in general there exist three critical points, 
which is illustrated in Fig. 2 for N=10 bosons for the limiting case of x — 
— y/7/2. The thick lines in Fig. 2 give points in the (CP) plane, which are 
local minima of the energy functional (2). The shaded area is the parameter 
range of (, where two local minima of the energy functional coexist. The lower 
dashed line gives the critical ( value where a deformed minimum appears, 
while the upper dashed line gives the critical point in ( where the spherical 
minimum disappears and only the deformed minimum is left. The dotted line 
gives the critical ( value where two coexisting minima are equally deep. The 
parameter region in between is small for any boson number. 
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Thus, as it is the aim of this work to identify the critical points in observables, 
and we do not expect to be able to distinguish between these three points (close 
lying in () in real nuclei, we restrict ourselves to the critical point given by 
condition (3) where the spherical minimum disappears, and which is given by 



The x dependence is just a finite N effect, and thus it is convenient to vary 
only the parameter ( for the investigation of phase transitions between vi- 
brational and rotational nuclei. Additionally we note that the choice of our 
parametrization has the convenient feature that in the large N limit we get 
Cc = 0.5. 

While, due to their physical meaning, the endpoints of the line of critical 
points between x — an d X — ~ V^/2 in Fig. 1 can be approximately related 
to the non-IBA symmetries E(5) and X(5), we see that a much richer structure 
shows up in the IBA, where critical points occur over the whole transitional 
region between these legs of the symmetry triangle. 

Since we are interested in obtaining signatures for critical points in observ- 
ables including matrix elements, we now survey the behavior of Q-invariants 
[18,19] in the transition regions. Recently, the concept of Q-invariants has 
been re-investigated in the framework of the IBA model and the Q-phonon 
approach [20,21], and the behavior of these moments across the gamut of nu- 
clear collective structures has been elucidated [22-24]. These invariants rep- 
resent quadratic and higher order moments of the quadrupole operator. The 
invariants are denoted q n and K n = q n /<h^ ■, an d are defined by expressions 
of the generic type 

Qn ~ (*o|Ql •Q2"-Qn|tfo) (5) 



where \l/o is the ground state wave function, and where intermediate angular 
momentum couplings in the operator are omitted for simplicity. 

For the IBA [8] , the Q-invariants have been evaluated over the entire symme- 
try triangle of Fig. 1. To show the extreme cases, we first focus on the two 
transition paths U(5)^SU(3) (x = ~V7/2) and U(5)<->0(6) (x = 0). We note 
that the invariants q 2 , K 3 , K4, and <r 7 = Kq — K\ represent, respectively, the 
quadrupole deformation, the triaxiality, the softness of the nuclear shape in 
(3, and in 7. 

We first study the U(5)^SU(3) transition and obtain the results shown for 
N=10 in the top row of Fig. 3 for q 2 , K 4 and cr 7 . Each of these exhibits a rapidly 
changing behavior which has a turning point ( t near ( = 0.5. To investigate this 
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in more detail, the second row of Fig. 3 shows the first derivatives with respect 
to (. Again there is a striking consistency of behavior: the first derivative has 
an extremum at essentially the same point for each invariant. 

Specifically, the turning points (the zeros of the second derivatives) are: ( t = 
0.54 for q2] Ct = 0-53 for K4, and ( t = 0.52 for <j 7 . In the coherent state 
formalism, for N=10, one obtains ( c = 0.54 for the U(5)^SU(3) case. This 
is very close to the turning points in q2,K± and <r 7 : that is ( t ~ Cc- This 
correspondence between the turning points and the critical points is the main 
result of this work. The small differences probably represent a finite boson 
number effect. 

This identification of a special point along the structural evolution from vi- 
brator to rotor is apparent even in the simplest observables as well. In Fig. 4 
we show the behavior of the structural observables R4/2 = E(Ai)/E(2f ) and 
B(E2 : 2f -> 0t) for the U(5)^SU(3) transition, again for N=10. Clearly, 
as seen in the first derivative plots in the second row, both quantities exhibit 
their steepest rates of change near the critical points. Here, the first derivative 
has an extremum at (t = 0.54 for both R4/2 and the B(E2) value. In this lat- 
ter case, this result is not surprising since this B(E2) value and q2 are directly 
related. 

The existence of three critical points on the U(5)^SU(3) transition path seems 
not to be reflected in the Q-invariants, which may be explained by the very 
compact parameter region in ( where these critical points occur, while the 
peaks in the derivatives have a certain width. Also note that fluctuations, 
resulting from the limited numerical accuracy of the Phint code used for 
these calculations, have been smoothed by the use of splines. Thus, perhaps 
the three critical points just cannot be resolved in the observables due to 
numerical truncations. 

Returning to the Q-invariants, similar results apply in the U(5)^0(6) region. 
Fig. 5 (left panels) shows this for q 2 and K4. In this case the turning points 
(determined from the rates of change), are: ( t = 0.60 for q 2 and ( t = 0.56 
for K 4 . From Eq. (2), the coherent state formalism gives ( c = 0.56 for N=10. 
Again the (t and ( c values obtained from the behavior of the Q-invariants and 
from the coherent state formalism are quite close. Lastly, we note that the 
rate of change of q2 and K4 in the U(5)^0(6) case is much less than in the 
first order U(5)— >SU(3) transition region. For example (dq2/d() ma x ~ 800 for 
U(5)^SU(3) while it is only -200 for U(5)<->0(6). Also, the widths of the first 
derivative curves are much wider (corresponding to a more gradual structural 
evolution) in the U(5)<-h>0(6) case. 

Using the IBA, it is also possible to investigate internal paths in the symmetry 
triangle. In particular internal straight line trajectories, starting from U(5), 
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will correspond to x values between and — \fl /2, allowing a full mapping of 
transitional trajectories. We illustrate such results by showing the change of 
the first derivative of the shape invariant K4 for various values of x i n Fig- 
6. The minima of the derivatives follow the line of critical points that is also 
given in the coherent state formalism, with only a small x dependence. 

Finally, in regard to Q-invariants, we look at the 0(6)<->-SU(3) transitional 
region. The right panels of Fig. 5 show the behavior of q<± and its derivative. 
Note that the shape is qualitatively different than in the other transition 
regions, showing a gradually asymptotic curve and a first derivative against x 
(the appropriate variable for this region) which is monotonic. No critical point 
is definable in this region of x values, except when 0(6) itself is reached (see 
ref. [11]). 

Another theme in nuclear structure recently has been the study of order and 
chaos for different structures. It was shown in ref. [25] that nuclear systems 
display ordered spectra at and near the three symmetry limits of the IBA, 
but that there is a rapid onset of chaotic behavior away from these bench- 
mark regions. (See Fig. 1 of ref. [25] but note that the symmetry triangle is 
differently defined therein.) Recently, Cejnar and Jolie [26,27] have developed 
the concept of wave function entropy as an alternate (and physically intuitive) 
way of studying the relative complexity of nuclear wave functions. Basically, 
the entropy of a state is a measure of its spreading within a given basis. Note 
that this is not the same as the chaoticity (which is basis invariant) since a 
wave function may have high entropy in one basis [e.g., U(5)] and low entropy 
in another [e.g., SU(3)]. 

Now that we showed a visible effect of critical points in various observables, it 
is interesting to see whether effects of a phase transition can also be seen in the 
wave functions and thus the wave function entropy. A rise of the wave function 
entropy can be expected in moving from one limit to another, but the question 
is whether it also appears in a close region with turning points which coincide 
with the turning points of the previously mentioned observables. Thus, we 
define [26] a quantity, called W®, for a state \1>, that can be written in the 
basis B as * = £™ a it3 \^ B >, as 

= K| 2 m|afJ 2 (6) 

ig=l 

where n is the number of basis vectors. If ^ coincides with a basis vector, then 
W® = 0. If ^ is uniformly spread out over the basis £>, then W® Inn. 

A physically intuitive expression of the entropy is the quantity [27] 

n?„ = exp W$ (7) 
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which expresses a kind of "effective number" of wave function components. 
For a "pure" state wf// = 1 and for a fully de-localized state nfjj ~ n. 

To properly normalize the entropies we define the entropy ratio 

b = ex P W E - 1 m 

exp(W GOE ) - 1 { J 



relative to that for the Gaussian Orthogonal Ensemble [27]. The ratio r% varies 
from for a pure (localized in the basis B) state to ~1 for a highly mixed 
state (see ref. [27] for a more detailed discussion of this normalization). 

We show the results in Fig. 7 for rj?+ and its derivative as a function of the order 

parameter ( for the U(5)^SU(3) and U(5)^0(6) transition regions (all for 
N=10). The entropy ratio for the ground state undergoes a very rapid change 
near ( c for both transition regions. We note that for larger boson numbers N 
the transition becomes much sharper (see Fig. 6 in [11]). For the U(5)^SU(3) 
and U(5)^0(6) phase transitions, it is easy to read the turning points, (t, 
values from the derivative plots, obtaining ( t = 0.52 and (t = 0.59 [in a U(5) 
basis], respectively, compared to values of ( c = 0.54 and ( c = 0.56 from the 
coherent state formalism. We note that the steepness of the entropy functions 
against ( increases with boson number N, as pointed out in ref. [28] . This also 
holds true for the observables studied above. 

To conclude, from the behavior of several rather different quantities, the Q- 
invariants, the simple observables R 4 / 2 and B(E2 : Z\ — > 0^), and the wave 
function entropy, we have shown that critical points of the phase transitional 
regions U(5)^SU(3) and U(5)^0(6) are reflected in the behavior of these 
observables along these evolutionary trajectories. This result was obtained 
for finite boson numbers, making it possible to investigate effects of valence 
particle number on the singularities. 

We are grateful to N.V. Zamfir, F. Iachello, J. Eberth and K. Heyde for useful 
discussions, and to P. Cejnar for the entropy calculations. Work supported 
by the U.S. DOE under Grant number DE-FG02-91ER40609 and by the DFG 
under Project number Br 799/10-1 and by NATO Research Grant no. 950668. 
One of us [RFC] is grateful to the Institut fur Kernphysik in Koln for support. 
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0(6) 
Ax. Asym. Rotor 
( y - soft) 




Fig. 1. Symmetry triangle of the IBA model. The U(5)<->0(6) leg is characterized by 
X = and varying (, while the U(5)^SU(3) transition region has \ = ~V7/2 and 
C is varied. The dashed line indicates the phase transitional region where critical 
points are found. 



0.56 




Fig. 2. The thick lines represent the locus in the ((,(3) parameter space where the 
energy functional of the coherent state formalism has a local minimum. The thick 
line at (3 = extends downwards to ( = 0. The results are shown for the case of 
N=10 bosons. Dashed lines mark critical ( values where one minimum disappears 
(the spherical one at and above the larger value, the deformed one at and below the 
lower value). Only in the shaded area two minima coexist. The dotted line marks 
the critical ( value where these two minima are equally deep. 
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U(5) -> SU(3) 
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Fig. 3. Behavior of (72, and <r 7 , and their first derivatives with respect to £, for 
the U(5)^SU(3) transition region, calculated for N=10 bosons. 



U(5) -> SU(3) 



50 




Fig. 4. Similar to Fig. 3 (for N=10) for the observables R 4/2 and B{E2 :2f^0f) 
for the U(5)^SU(3) transition region. 
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Fig. 5. Similar to Fig. 3 (for N=10), for qi and K4, for the U(5)<->0(6) transition 
region (left panels), and for qi in the 0(6)<-^SU(3) (note here with respect to x) 
transition region (right). 




Fig. 6. First derivative of K4 (for N=10), but for various values of the parameter %. 
A peak indicating a phase transition occurs for every value of %. The dependence 
of its position on x is a finite N effect. 
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U(5) -> SU(3) 



U(5) -> 0(6) 




Fig. 7. The entropy ratio (for N=10) for the 0^ state (top row) in the three transition 
regions, plotted against £ and given, for each region, in two bases as indicated [e.g., 
U(5) and SU(3) for the U(5)^SU(3) transition]. The lower panels give the derivative 
of the entropy ratio against £ in the appropriate basis. 
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